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Abstract 

With a new detector setup and the high-resolution performance of the fragment separator FRS at GSI we 
discovered 57 new isotopes in the atomic number range of 60< Z < 78: 159_161 Nb, 16 °- 163 p m; 163_166 Sm, 

167-168 EUj 167-171 Gd; 169-171^ 171-174^ 173~176 HOj 176-178^ 178-181^ 183-185 Ybj 187-188 LUj 191 Hf? 
193-194 Taj 196-197 Wj 190-300^ 201-203 OSj 204-205j r &nd 206-209p t Xhe new isotopes have been unam _ 

biguously identified in reactions with a 238 U beam impinging on a Be target at 1 GeV/u. The isotopic 
production cross-section for the new isotopes have been measured and compared with predictions of differ- 
ent model calculations. In general, the ABRABLA and COFRA models agree better than a factor of two 
with the new data, whereas the semiempirical EPAX model deviates much more. Projectile fragmentation is 
the dominant reaction creating the new isotopes, whereas fission contributes significantly only up to about 
the element holmium. 

Key words: NUCLEAR REACTIONS Be( 238 U,A), E = 1 GeV/u; measured Nb-Pt fragments isotopic 
production cr; Comparison with previous results and model predictions. 
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1. Introduction 

Heavy neutron-rich nuclides are of great interest for 
nuclear astrophysics and basic nuclear spectroscopy. 
This becomes immediately obvious when one looks 
at the predicted path for r-process nuclei and their 
decay. The study of shell evolution far off stability 
and towards the expected magic numbers iV=82, 126 
and thus the waiting points of the r-process nuclides 
are of interest for both fields. The accurate knowl- 
edge of the atomic masses and lifetimes are essential 
for the understanding of the nucleosynthesis 0, 0, HJ . 
Presently, the corresponding theories when applied 
to newly opened experimental territories still deviate 
significantly from the results of measurements. 

Experimentally the area of heavy neutron-rich nu- 
clides is difficult to reach because of the low produc- 
tion cross-sections and the great challenge of sepa- 
ration and isotopic identification. Relativistic ener- 
gies of the reaction products and the high ion-optical 
resolution of the in-flight separator FRS |4| are the 
keys to the frontiers in this domain of nuclides. High 
velocities are required to reduce the number of popu- 
lated ionic charge-states for each element, mainly to 
bare fragments with a low contamination of H- and 
He-like ions. 

Several milestones in nuclear physics have been 
achieved with the FRS, like the discovery and spec- 
troscopy of the double magic nuclei 100 Sn 0, Q or 
the discovery of 2p radioactivity in 45 Fe In the 
pioneering experiments with uranium projectile fis- 
sion with the FRS at 750 MeV/u more than 120 
new isotopes have been discovered, among them 78 Ni 
8]. These achievements have launched a new re- 
search activity for fission studies and an area of 
applied physics towards accelerator-driven reactors 
and nuclear- waste transmutation (lp| . Owing to the 
success of in-flight fission at high energies, all next- 
generation in-flight exotic nuclear beam facilities in- 
clude the production via projectile fission. For ex- 
ample, very recently the new powerful Radioactive 
Ion Beam Factory (RIBF) in RIKEN has success- 



fully started the experimental program with the dis- 
covery of new isotopes in the atomic number range 
of 26 < Z < 56 produced via in-flight fission of an 
intense 238 U beam ll|. The search for the neu- 
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tron dripline at low Z has been a major research 
activity at GANIL (France) [H, NSCL (USA)[H 
and RIKEN (Japan) 0, [l2| for many years. At the 
NSCL facility at MSU (USA) new neutron-rich iso- 
topes have been observed via reactions with a 48 Ca 
and 76 Gc beam at about 140 MeV/u 0[l3|. 

In the recent years, the intensity for 238 U beams 
provided by the GSI accelerators have increased al- 
most by a factor 10, which has opened new perspec- 
tives for the production and study of the heaviest 
projectile fragments 16|, [l7[ . Even along with mass 
measurements at the FRS-ESR facility new isotopes 
have been observed [il[l7|. Thefra gmentation reac- 
tion of 208 Pb seems also to be very promising for pro- 
duction of neutron-rich isotopes as it has been proven 
in the recent FRS experiments 18, 19j . 

In this letter we report on the discovery of 57 new 
neutron-rich isotopes in the element range of Nd to 
Pt at the FRS by the use of new particle identification 
methods [Ej]. 

2. Experimental technique 

The experiment was performed with the SIS-18 
synchrotron of GSI Darmstadt, which delivered a 1 
GeV/u 238 U beam in spills lasting 0.5-2 s with a rep- 
etition period of 2-4 s. The beam impinged on a 
1.6 g/cm 2 thick beryllium target placed at the en- 
trance of the projectile FRragment Separator (FRS) 
4] . The primary beam intensity was of the order of 
2 x 10 9 ions/spill. The 238 U intensity was recorded 
by a calibrated secondary-electron transmission mon- 
itor [22|. The reaction products were separated by 
the FRS operated in an overall achromatic ion-optical 
mode. A schematic view of the FRS and the ex- 
perimental setup is shown in Fig. [T] The spatial 
separation in flight was achieved by twofold appli- 
cation of the Bp-AE-Bp method, i.e., the atomic 
energy losses in two degraders, located at the first 
(Fl) and second (F2) focal planes, were measured 
via magnetic rigidity analysis. In this way, the re- 
action products are spatially separated and by the 
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use of various detectors their nuclear charge Z and 
mass number A could be determined. After the first 
magnetic selection and the 2.5 g/cm 2 thick aluminum 
degrader at Fl, the reaction products were slowed 
down in an aluminium disk degrader located at the 
intermediate focal plane F2. With the disk angle the 
degrader shape was tuned to preserve the achroma- 
tism. Even at these relativistic velocities the atomic 
charge states of the heavy fragments represent a sep- 
aration problem. Therefore, medium Z-material nio- 
bium foils were placed both behind the target and 
the F2-degrader to enhance the yield of bare frag- 
ments. The thicknesses of these electron strippers 
were 223 mg/cm 2 for the first and 106 mg/cm 2 for 
the latter. The total thickness of the F2-materials 
including detectors was 1.43 g/cm 2 aluminium equiv- 
alent. 

The complete particle identification in-flight was 
performed on an event-by-event basis with time-of- 
flight (ToF), energy-deposition (AE 1 ) and magnetic 
rigidity measurements (ToF- AE' -Bp method). The 
ToF measurement was performed with two plastic 
scintillator detectors, one located at F2 and the other 
one at the final focal plane (F4). The flight path was 
about 37 m between the two detectors. The ToF 
value for the selected isotopes was of the order of 160 
ns in the laboratory frame. At the exit of the FRS, 
two ionization chambers (MUSIC) filled with P10 gas 
at 1 atm. pressure were mounted with a 104 mg/cm 2 
copper stripper placed in between. The MUSIC de- 
tectors [23| delivered the energy-deposition signals of 
fragments, thus providing the information of their 
atomic numbers. The velocity dependence of the 
penetrating ions was taken into account before the 
energy-deposition signal was applied for Z identifica- 
tion. The magnetic rigidity measurements were per- 
formed with four time-projection chambers (TPC), 
two located at the dispersive focal plane (F2) and two 
others mounted at the exit of the FRS. The TPC pro- 
vided full tracking information (angle and position) 
for the transmitted fragments. The event- by- event 
identification, and thus also the in-flight separation, 
were verified by the isomer tagging technique 21 1. 



incidence with the incoming ions. With these mul- 
tiple redundant measurements and the two-stage in- 
flight separation criteria, we achieved an unambigu- 
ous isotope identification. Finally, a range selection 
was applied in addition because the selected frag- 
ments were stopped in a layer of matter viewed ei- 
ther by the RISING germanium detector setup 24(, 
composed of fifteen Euroball cluster of seven crystals 
used in Air configuration, or the simpler isomer tag- 
ging device [21(, which consisted of the two electro- 
mechanically cooled Ge detectors, a stopper foil, and 
a veto scintillation counter. 
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In the range of the particle identification spectrum, 
known fis isomers were selected ( 171m Tm, 172m Yb 
and 175m Lu) whose gamma rays were recorded in co- 



Figure 1: (color online) Schematic view of the FRS and the ex- 
perimental setup. The magnetic dipole and quadrupole mag- 
nets of the FRS and the target and different focal-plane ar- 
eas are depicted. The different energy degraders at the focal 
planes Fl, F2, and F4 are schematically shown. The detector 
setups placed in the central and final focal planes are shown 
as zoomed inserts in the picture. After identification in flight 
the ions are stopped in an implantation detector viewed by a 
gamma detector array. 



3. Data analysis 

In the experiment four different Bp settings of the 
FRS were applied, which were chosen to yield op- 
timum intensities for bare 172 Dy, 194 0s, 198 0s and 
202 Os ions, i.e., the magnetic fields were selected to 
center these isotopes during their travel inside the 
FRS system. The data collected in each setting were 
processed by the following procedure based on a com- 
bination of ToF, position and energy-deposition. 

In the first step of analysis a two-dimensional plot 
of the energy deposition in the first MUSIC detector 
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as a function of the energy deposition signal from the 
second MUSIC was created. Only bare ions in both 
detectors were chosen for further analysis, thus re- 
jecting ions with different charge states in both MU- 
SICs or secondary reaction products created in the 
detector material. 

In the second step, a distribution of the recon- 
structed nuclear charge as a function of the energy 
loss of the ions in the degrader located at F2 was cre- 
ated. This allowed a clear identification of the group 
of ions which did not change their charge state while 
penetrating the matter placed at F2. 

Unambiguous isotope identification requires an ad- 
ditional selection based on the ions' position at the 
final focus (Bp), energy deposition signals in both 
scintillators and a correlation between the measured 
angles at the intermediate and final focal planes. The 
latter condition reflects the ion-optical image condi- 
tions which are needed also for the correct Bp deter- 
mination. An example of the final isotope identifica- 
tion plot is shown in Fig. [3J The projection of this 
plot on the A/q axis selecting every element covered 
in the Bp settings is given in Fig. [3] 

4. Production cross-sections 

The data recorded for each Bp setting were ana- 
lyzed by the described procedure to achieve an unam- 
biguous isotope identification. In the next step the 
production cross-sections of individual isotopes were 
determined according to: 



N 



f 



T op t T sec Pa YNp fdt 



(1) 



where Nt is the number of registered ions of a certain 
isotope, T opt the ion-optical transmission, T sec a cor- 
rection for secondary reactions in the matter placed 
after the target (e.g., degraders and detectors), Pq the 
probability that an ion remains fully stripped in both 
stages of the separator, Y the correction for losses of 
primary beam and fragments due to nuclear reactions 
in the target material, N p the total number of 238 U 
ions and fdt the correction for the dead-time losses 
of the data acquisition system. All secondary reac- 
tions in the target and the matter in the focal planes 
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Figure 2: (color online) Identified atomic number Z of the in- 
coming ions as a function of their A/q ratio, at the final focal 
plane (F4). In the A/q ratio, the time-of-flight and magnetic 
rigidity analysis information is included. The plot shows the 
superimposed data acquired in all Bp settings. The solid line 
shows the border of hitherto unobserved isotopes, i.e., the dis- 
covered new isotopes are on the right hand side of the border 
line. The black rectangles show the group of ions correspond- 
ing to the Bp setting of the spectrometer. 



are taken into account by applying the Benesh-Cook- 
Vary formula [25| . 

Calculations of charge-state distributions were per- 
formed for the heavy fragments of interest with the 
CHARGE code [26( yielding P values in the range 
of 0.9 and 0.6 for Nd-Pt isotopes, respectively. The 
ion-optical transmissions T opt have been calculated 
by using the Monte-Carlo simulation program MO- 
CADI [27|, l28( assuming the kinematics of projectile 
fragmentation [35| . The transmission values were ob- 
tained separately for each Bp setting. The typical 
T op t values were of the order of 0.4-0.6 for isotopes 



with A/q close to the reference setting of each FRS 
setting. This relative small value is due to the tight 
slits settings of the FRS which were applied in order 
to decrease the number of contaminants, mainly light 
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Figure 3: (color online) Projection of the identification plot shown on Fig.[2]constructed for all elements covered by the different 
Bp settings of the FRS. For each element the arrow indicates the lightest of the isotopes (marked by the mass numbers) observed 
for the first time. 
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fission fragments (Z<60), reaching the F4 area. 
5. Results 

In this experiment we discovered 57 new isotopes 
with atomic numbers in the range 60 < Z < 78: 

159-161]^ 160-163p m 163-166g m 167-168g u 
167-m Gd ' 169-171 Tb ' 171-174 Dy) ' 173-176 Ho ' 
176-178 Erj 178-181 Tm) 183-185 Yb; 187-188^ 191^ 
193-194rp a 196-197^ 199 - 200 Re 201-203Q g 204-205j r 

and 206 - 209 Pt and measured their production cross 
section. The new isotopes have been unambiguously 
identified in flight by applying a two-fold Bp-AE-Bp 
separation scenario and redundant ToF-AE'-Bp 
analysis using new detector systems [21(. The 
observed new isotopes are presented for the different 
elements in Fig. [3] For each element the arrow 
indicates the lightest of the isotopes observed for 
the first time. For all isotopes discovered in this 
work we can set a lower limit of half-life to 300 ns 
which corresponds to the time of flight between the 
production target and F4. 

The isotopic production cross-section have been de- 
termined in this experiment down to the pico-barn 
level. Particularly, the steep descent of the yields 
of the neutron-rich isotopes has been mapped. The 
cross-section obtained for the production of Nd-Pt 
isotopes are shown in Fig. [4] and listed in Table [TJ 
The error bars are mainly determined by the statis- 
tics due to tiny production cross-sections, partially 
due to the low transmission of isotopes character- 
ized with a Bp value far from the central field setting 
with respect to the optical axis. The latter aspect 
has been partially taken care of by the different field 
settings applied in these measurements. The mea- 
sured cross-sections in this experiment represent the 
outskirts of the chart of nuclides. Therefore, it is ob- 
vious that the present results cannot be compared 
with previous experimental data. However, a good 
orientation for the consistencies of the cross-section 
evolution can be obtained with the comparison to 
results from previous experimental studies [sT 



31|. 



Although the latter data from the literature are ob- 
tained from reactions of 1 GeV/u 238 U beam with 
a liquid hydrogen target the continuous transition of 
the two experimental data sets is remarkable. The 



excitation in the reaction with hydrogen target nu- 
clei should be lower than in our case with a beryl- 
lium target but this plays obviously a minor role for 
the compared fragment distribution in the overlap 
region. Our cross-section are compared with calcula- 
tions based on the ABRABLA [II] and COFRA [H] 
models and shown in Fig. |4j The ABRABLA model 
is an improved abrasion-ablation model which takes 
also into account microscopic structural effects and 
the contribution of fission. COFRA and EPAX have 
both the fission process not included. The predictions 
of the semiempirical EPAX parametrization 34] arc 
given as well. 

Since ABRABLA is a Monte-Carlo code, a long 
computation time is required to reach very-low pro- 
duction cross-section. Thus COFRA, the analyti- 
cal version of the abrasion-ablation model, was used 
for the most neutron-rich nuclei. The predictions of 
ABRABLA show an overall good agreement with the 
experimental data. The used version of the code in- 
volves fission as a possible deexcitation step of the 
created prefragment. The lighter isotopes in the el- 
ement range Nb-Dy show a significant contribution 
from the fission process which starts to play a minor 
role for heavier nuclei investigated in this work. The 
latter statement was deduced from the ABRABLA 
calculation and is also visible by the underestimation 
of EPAX, which takes into account only projectile 
fragmentation. 

The COFRA model agrees well with the experi- 
mental data, however, a small gradual deviation be- 
tween results of the measurements and the model pre- 
dictions can be observed for the heaviest (Os-Pt) nu- 
clei. For most of the cases the overestimation is less 
than a factor of two compared with the measured 
values which means it is well suited to make valid 
predictions for new studies in this field of neutron- 
rich isotopes. The EPAX predictions largely overes- 
timates the production cross-sections in the region of 
neutron-rich nuclei in the element range of Yb to Pt. 
This reflects the present limitation of that semiem- 
pirical model with its parametrization based on pre- 
vious data, hence our results could be a base for a 
new parametrization of EPAX. 
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Figure 4: (color online) Measured production cross-sections of fragments produced in the reaction 238 U (1 GeV/u) + Be (black 
circles), shown together with the experimental results obtained by Bernas et al. l2f| (Z = 60 — 64), [3fi (Z = 65 — 73) and 
Taieb et al. [H[ (Z = 74 - 78) in the reaction 238 U (1 GeV/u) + p (blue open circles). The red dashed line represents the 
predictions of the ABRABLA model [13 and the continuous red line shows the results of COFRA (Z = 73 - 78). The 
green dotted line shows the prediction of EPAX model [34H . For each element the arrow indicates the lightest of the isotopes 
observed for the first time. 
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Tabic 1: Production cross-section of new isotopes measured in this work in reaction 1 GeV/u 238 U on Be. 



Isotope 


cr(nb) 


Isotope 


cr(nb) 


Isotope 


ff(nb) 


Isotope 


cr(nb) 


lb9 Nb 


20(6) 


itiy Gd 


3.2(5) 


lW Er 


2.7(0.3) 




0.0034(17) 


160 Nb 


3.9(1.4) 


170 Gd 


1.9(8) 


178 Er 


1.1(2) 


199 Re 


0.0076(27) 


161 Nb 


0.7(4) 


171 Gd 


0.6(4) 


178 Xm 


10.6(1.8) 


200 Re 


0.0025(14) 


160p m 


180(15) 


i69 Xb 


71(4) 


179 Xm 


2.1(3) 


201 Os 


0.033(7) 


161p m 


45(6) 


170 Xb 


14.6(1.1) 


180 Xm 


0.79(17) 


202 Os 


0.0044(2) 


162p m 


17(6) 


171 Xb 


3.3(0.4) 


181 Tm 


0.11(6) 


203 Os 


0.0035(17) 


163p m 


2.2(1.0) 


171 Dy 


49(3) 


183 Yb 


0.21(6) 


204 T r 


0.017(5) 


163 Sm 


53(4) 


172 Dy 


15.0(1.1) 


184 Yb 


0.028(9) 


205 j r 


0.0026(15) 


164 Sm 


14(2) 


173 Dy 


3.5(4) 


185 Yb 


0.007(3) 


206p t 


0.033(11) 


165 Sm 


8.3(4) 


174 Dy 


0.57(18) 


187 Lu 


0.043(9) 


207p t 


0.008(3) 


166 Sm 


0.7(5) 


173 Ho 


65(6) 


188 Lu 


0.010(3) 


208p t 


0.0027(15) 


167 Eu 


2.9(7) 


174 Ho 


12(1) 


191 Hf 


0.0043(25) 


209p t 


0.0020(14) 


168 Eu 


0.7(3) 


175 Ho 


3.4(4) 


193 Ta 


0.017(5) 






167 Gd 


71(4) 


176 Ho 


0.51(14) 


194 Ta 


0.0037(19) 






168 Gd 


14.9(1.2) 


176 Er 


1.19(12) 


196 W 


0.018(4) 







6. Summary 

The results of the present experiment opens a 
new field for nuclear spectroscopy and also for nu- 
clear astrophysics in the heavy nuclei range. With 
a new detector setup and the high-resolution per- 
formance of the fragment separator FRS we discov- 
ered 57 new isotopes in the atomic number range of 
60 < Z < 79. The new isotopes have been unambigu- 
ously identified in projectile fragmentation reactions 
with a 238 TJ beam at 1 GeV/u. The isotopic pro- 
duction cross-section for the new isotopes have been 
measured and compared with the predictions of the 
ABRABLA, COFRA, and EPAX models. In gen- 
eral, with COFRA a good overall agreement has been 
achieved to be confident for reliable predictions in un- 
known territory. The next steps in this experimental 
campaign will be half-life and mass measurements, 
as well as decay spectroscopy after implantation in 
silicon detectors. 
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